Applications of quark-hadron duality in F2 structure function 
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Inclusive electron-proton and electron-deuteron inelastic cross sections have been measured at 
Jefferson Lab (JLab) in the resonance region, at large Bjorken x, up to 0.92, and four-momentum 
transfer squared Q 2 up to 7.5 GeV 2 in the experiment E00-116. These measurements are used 
to extend to larger x and Q 2 precision, quantitative, studies of the phenomenon of quark-hadron 
duality. Our analysis confirms, both globally and locally, the apparent "violation" of quark-hadron 
duality previously observed at a Q 2 of 3.5 GeV 2 when resonance data are compared to structure 
function data created from CTEQ6M and MRST2004 parton distribution functions (PDFs). More 
importantly, our new data show that this discrepancy saturates by Q 2 ~ 4 Gev 2 , becoming Q 2 
independent. This suggests only small violations of Q 2 evolution by contributions from the higher- 
twist terms in the resonance region which is confirmed by our comparisons to ALEKHIN and 
ALLM97. We conclude that the unconstrained strength of the CTEQ6M and MRST2004 PDFs 
at large x is the major source of the disagreement between data and these parameterizations in the 
kinematic regime we study and that, in view of quark-hadron duality, properly averaged resonance 
region data could be used in global QCD fits to reduce PDF uncertainties at large x. 

PACS numbers: 25.30.Fj, 24.85.+p 



I. INTRODUCTION 



'Deceased To understand how Quantum ChromoDynamics 

(QCD) works remains one of the great challenges in nu- 
clear physics today. The challenge arises from the fact 
that the degrees of freedom observed in nature, hadrons 
and nuclei, are not the same as the ones appearing in the 
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QCD Lagrangian, quarks and gluons. The challenge is 
then to formulate a connection between the description 
of hard, or short-distance, scattering processes which can 
be calculated perturbatively in terms of quark and gluon 
degrees of freedom and their weak couplings, and soft, 
or long-distance, scattering processes, where the physi- 
cal asymptotic states are prominent and the quarks and 
gluons interact strongly. 

Given these strong quark-gluon interactions, or the 
large value of the strong coupling constant a s , the spec- 
tra of the asymptotic hadron states are not calculable 
within a perturbative QCD (pQCD) framework, and are 
difficult to directly connect to the underlying quark- 
gluon or parton dynamics. Yet, several instances ex- 
ist in nature where the behavior of low-energy scatter- 
ing cross sections, averaged over appropriate energy in- 
tervals, closely coincide with asymptotically high-energy 
scattering cross sections, calculated in terms of quark- 
gluon degrees of freedom. This phenomenon is referred 
to as quark-hadron duality, and may be a general prop- 
erty of quantum field theories with inherent weak and 
strong coupling limits, with QCD as a prime example. 

The observation of a non-trivial relationship between 
inclusive electron-nucleon scattering cross sections at low 
energy, in the region dominated by the nuclcon reso- 
nances, and that in the deep inelastic scaling regime at 
high energy predates QCD itself. While analyzing the 
data from the early deep inelastic scattering experiments 
at SLAC, Bloom and Gilman observed [3, HJ that the 
inclusive structure function at low hadronic final state 
mass, W, generally follows a global scaling curve which 
describes high-PF data, and to which the resonance struc- 
ture function averages. Following the development of 
QCD in the early 1970s, Bloom-Gilman duality was re- 
formulated in terms of an operator product (twist) ex- 
pansion of moments of the structure functions [H, 0|. 
This allowed a systematic classification of terms respon- 
sible for duality and its violations in terms of so-called 
higher-twist operators which describe long-range inter- 
actions between quarks and gluons. However, this de- 
scription could not explain why particular multi-parton 
correlations were suppressed, and how the physics of res- 
onances gave way to scaling [5[ . 

Since then, with the development of high luminosity 
beams at modern accelerator facilities such as JLab, a 
wealth of new information on structure functions, with 
unprecedented accuracy and over a wide range of kine- 
matics, has become available. One of the striking findings 
of the new JLab data [ESQ * s * na * Bloom-Gilman dual- 
ity appears to work exceedingly well, down to Q 2 values 
as low as 1 GeV 2 or even below. This is considerably 
lower than previously believed, and well into the region 
where a s is relatively large. Furthermore, the equivalence 
of the averaged resonance and scaling structure functions 
appeared to hold for each resonance, over restricted re- 
gions in W, so that the resonance-scaling duality holds 
also locally. It was also found that quark-hadron duality 
manifests itself in the separated proton transverse (F[) 



and longitudinal (F^) structure functions. 

The more recent JLab resonance structure function 
studies have revealed an important application of dual- 
ity: if the workings of the resonance-deep inelastic inter- 
play are sufficiently well understood, the region of high 
Bjorken-x [x > 0.7, where x is the longitudinal momen- 
tum fraction of the hadron carried by the parton in the 
infinite momentum frame) would become accessible to 
quantitative studies. This region remains largely unex- 
plored experimentally due to the requirement of high- 
energy beams with sufficiently high luminosity. 

The x — ► 1 region is an important testing ground 
for nonperturbative and perturbative mechanisms un- 
derpinning valence quark dynamics, and is vital to map 
out if we hope to achieve a complete description of nu- 
cleon structure. Data from the nucleon resonance region, 
where quark-hadron duality has been established, could 
be used to better constrain QCD parameterizations of 
parton distribution functions (PDFs), from which also 
the hadronic backgrounds in high-energy collisions are 
computed The large- a; region also constitutes an ap- 
preciable amount of the moments of polarized and un- 
polarized structure functions, especially for the higher 
moments. It is precisely these moments that can be cal- 
culated from first principles in QCD on the lattice [9j, in 
terms of matrix elements of local operators. 

Note that, since the x dependence of structure func- 
tions cannot be calculated on the lattice directly, one 
cannot easily use the lattice to learn about the degree to 
which duality holds locally. Indeed, the ability to calcu- 
late a leading-twist moment on the lattice implicitly uses 
quark-hadron duality to average the resonance contribu- 
tions to a smooth, scaling function. 

In this article, we quantitatively study the application 
of quark-hadron duality to access parton dynamics in the 
region of large x, up to x ~ 0.9. For this, we accumu- 
lated a series of inclusive electron-proton and electron- 
deuteron scattering data in the nucleon resonance region 
(W 2 < 4 GeV 2 ), at the highest momentum transfers ac- 
cessible at JLab. These data are at values of Q 2 far 
above where duality was quantitatively found to be valid 
in previous JLab experiments. The extracted Fi struc- 
ture function data are also compared with various state- 
of-the-art parameterizations of Fi world data to improve 
our understanding of parton dynamics at large values of 
Bjorken x. 

The article is structured in five sections. Section 2 
summarizes techniques of modeling the dynamics of the 
nucleon in terms of structure functions computed from 
PDFs and examines in detail few representative parame- 
terizations of the nucleon Fi structure function focusing 
on the large x region. Section 3 is an overview of the 
experimental apparatus utilized to collect these experi- 
mental data and of the analysis steps taken to extract the 
cross section and the F2 structure function. In Section 4 
we present our studies of the application of quark-hadron 
duality to gain insight in the parton dynamics at large x. 
In Section 5 we draw conclusions. 
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II. F 2 PARAMETERIZATIONS AT LARGE 
BJORKEN x 



The purpose of this Section is to give an overview of the 
techniques typically employed to map out the dynamics 
of the nucleon via structure functions. This discussion 
points out the importance, but also the difficulty, of ob- 
taining a parameterization of the F% structure function 
for the entire kinematic range. In particular, the ex- 
clusion of data in regions where the perturbative QCD 
mechanisms are not the only ones expected to contribute, 
greatly limits the applicability of these parameterizations 
and also our knowledge of the nucleon structure. In this 
context, quark-hadron duality might be the tool which 
could open kinematic regions not easily accessible oth- 
erwise to detailed studies. Four representative param- 
eterizations will be examined in detail with an empha- 
sis on the large x region: ALLM97 [nj, CTEQ6M [HI, 
MRST2004 nj and ALEKHIN EJ 



14[. These parame- 



terizations were used in our duality studies which will be 
presented in Sect. 4. 

Lastly, the parameterization of the structure function 
F p from Bourrely et al. [IH will also be considered. 
This parameterization is obtained from parton distribu- 
tion functions constructed in a statistical physical pic- 
ture of the nucleon, where the nucleon is viewed as a gas 
of massless partons (quarks, antiquarks and gluons) in 
equilibrium at a given temperature in a finite size vol- 
ume. The x dependence of the parton distributions is 
chosen to correspond to a Fermi-Dirac distribution for 
quarks and antiquarks, and to a Bose-Einstein distribu- 
tion for gluons. The parameterization involves a total 
of eight free parameters which are constrained by fitting 
high W 2 data from various experiments: NMC, BCDMS, 
E665, ZEUS, CCFR. A comparison of F p from Bourrely 
et al. to results from CTEQ6M and ALLM97 will be 
shown in Sect. 4. 

Empirical Parameterization of F p : ALLM97. 
ALLM97, proposed as an update of ALLM [Tj| published 
in 1991, is a Regge motivated parameterization extended 
to the large Q 2 regime in a way compatible with QCD 
expectations. The data set used to obtain the ALLM97 
fit coefficients included all 7*p measurements published 
up to 1997, with W 2 > 3 GeV 2 , and covering a wide 
range in Q 2 , < Q 2 < 5000 GeV 2 . The ALLM97 fit 
function has a total of 23 parameters, half of which are 
needed for the description of the low W 2 (high x) region 
where higher- twist terms are expected to be important. 
There are two important aspects to be noted in relation 
to the behavior of the ALLM97 parameterization at large 
x (see Figs. 1 and 2). On one hand, the data set used to 
obtain the fit coefficients is selected with a rather low W 2 
cut. Thus, it is expected that, if duality holds globally, 
the extrapolation of ALLM97 below W 2 of 3 GeV 2 in the 
resonance region will work reasonably well, on average. 
On the other hand, ALLM97 being an empirical fit, some 
of its shortcomings like unconstrained x and Q 2 depen- 
dence or inability to fully account for target mass effects 
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FIG. 1: (Color online) Existing F$ parameterizations at a Q 2 
value of 4 GeV 2 . CTEQ6M [tU, MRST2004 [TJ], ALLM97 
and ALEKHIN [HHi 

were used for the quark-hadron 
duality studies presented in Sect. 4. Target mass corrections 
were added to CTEQ6M and MRST2004 (see text). The M.E. 
Christy [171 ] parameterization is shown for comparison. 



will become obvious as we probe kinematic regimes out- 
side its domain of applicability. This will most likely be 
revealed in a clear manner when extrapolating to low W 2 
regions. 

QCD Parameterization of F p . Starting from two 
basic ideas of pQCD, factorization and evolution, the 
F% structure function can be calculated from PDFs ex- 
tracted from hard-scattering data [HI]. The theorem of 
factorization of long-distance from short-distance depen- 
dence in deep inelastic scattering (DIS) allows for the 
structure function to be expressed as a generalization of 
the parton model results: 



F2 P (x,Q 2 ) 



E 

=f,f,G 



Jo 4 

(1) 

where i denotes a sum over all partons (quarks, anti- 
quarks and gluons) inside the proton, Cj 1 are coefficient 
functions independent of the long-distance effects while 
<fii/p are parton distributions sensitive to the nonpertur- 
bative, long-distance effects inside the proton [l9| . The 
evolution, on the other hand, enables the systematic, per- 
turbative computation of logarithmic scale-braking ef- 
fects and ensures that measuring F2 P (x,Q 2 ) is enough 
to predict not only F^ p (x,Q 2 ) but also F^ p (x,Q 2 ) for 
all Q 2 , assuming that both Q 2 and Q 2 are large enough 
that a perturbative expansion in a s is still appropriate. 
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FIG. 2: (Color online) Existing F^ parameterizations at a Q 2 
value of 8 GeV 2 . CTEQ6M [n]], MRST2004 p3|, ALLM97 
[l(J and ALEKHIN [II [3 were used for the quark-hadron 
duality studies presented in Sect. 4. Target mass corrections 
were added to CTEQ6M and MRST2004 (see text). The M.E. 
Christy [T| parameterization is shown for comparison. 

This is typically done by using the DGLAP evolution 
equations to evolve the parton distributions: 

E / 1 f^(f,« s (Q 2 ))x^ /p (e,Q 2 ), (2) 

where Py are the evolution kernels (splitting functions) 
given by a perturbative expansion in a s , beginning with 
the leading order (LO) 0(a s ) but also calculable to 
higher orders, next-to- leading order (NLO) or next-to- 
next-leading order (NNLO). The kernels have the physi- 
cal interpretation as probability densities of obtaining a 
parton of type i from one of type j carrying a fraction of 
the parent parton's momentum. 

Thus, the three basic quantities are the coefficient 
functions C^ l \ the evolution kernels Py and the PDFs 
4>i/p. Of these, the first two are computed perturbatively 
as power series in a s . The physical non-perturbative par- 
ton distributions are extracted by combining theory and 
experiment and performing QCD fits. In a typical QCD 
fitting procedure the x dependence of the parton distri- 
butions is parameterized at some low scale, Qjj, where 
higher-order corrections in a s are expected to be negligi- 
ble, and then a fixed order ( either LO or NLO or NNLO) 
DGLAP evolution is performed to specify the distribu- 



tions at higher scales where data exist. A global fit to the 
data then determines the parameters of the input distri- 
butions. There is considerable freedom in choosing the 
parametric form of the parton distributions at scale Q 2 , 
|18j |. The parameterization should be general enough to 
accommodate all possible x. A typical choice is: 

<t>{x,Ql)=A Q x A ^l-x) A *P{x), (3) 

where P{x) is a smooth function of x, x Al and (1 — x) A2 
determines the small and large x behavior, respectively, 
and ^4o,i,2 are coefficients to be determined from fits to 
data. 

When performing QCD fits, there are several concep- 
tual difficulties to take into account. First of all, a QCD 
analysis of F2 measurements involves the use of the gluon 
distribution which is, a priori, unknown. In fact, the 
gluon distribution is the most uncertain of the PDFs and 
is particularly ill-determined for x > 0.3, with uncertain- 
ties reaching 200% by x = 0.5 [ll|. This in turn trans- 
lates in an uncertainty in the a s determination from QCD 
analysis of PDFs since there is a correlation between the 
hardness of the gluon and the magnitude of Aqcd, the 
quantity which sets the scale for ct s [2(j. The xF$ mea- 
surements should be able to provide a precise value of 
h-QCD since the gluon distribution does not enter into 
the evolution. However, the experimental uncertainties 
of xF$ are still larger than those of F2 and discrepan- 
cies in extracting xF$ were observed between different 
experiments in the region of x > 0.4 (2lj . 

Then, it should be pointed out that, a standard QCD 
analysis of PDFs does not take into account all residual 
Q 2 effects arising, for example, from higher-order radia- 
tive corrections in a s or/and non-perturbative higher- 
twist corrections. In particular, the higher- twist terms 
are non-factorisable and process dependent and QCD has 
no rigorous prescription to account for it. As a result, 
most PDFs extractions are performed using safe kine- 
matic cuts for data selection in order to exclude regions 
where higher- twist or/and higher-order corrections in a s 
play an important role. A typical set of cuts employed 
to select data for fitting is Q 2 > 2 GeV 2 and W 2 > 12 
GeV 2 (this type of W 2 cut rejects the whole resonance 
region). Limiting the data coverage to a particular range 
in Q 2 and W 2 will result in a limitation in the x cov- 
erage. For example, for a fixed Q 2 , a W 2 cut of W 2 > 
w u m wil1 limit the x ran S e to x < Q 2 /{Wf im - M 2 + 
Q 2 ). Considering that the x dependence of the PDFs is 
parameterized empirically, as exemplified in Eq. 3, and 
that the parameterization coefficients are extracted from 
fits to data, these data selection cuts, though they make 
possible the extraction of PDFs without the complica- 
tions specified above, yield to unconstrained strengths of 
the PDFs at large x [n| . 

To date, the large experimental uncertainties in the 
large x regime, when excluding the low W 2 data, pre- 
vent answers to basic questions as to why the d quark 
distribution appears to be softer than the u quark distri- 
bution. For the same reason, the d/u behavior at large 
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x, a critical test of the mechanism of spin-flavor sym- 
metry breaking, is highly unconstrained. Furthermore, 
knowledge of PDFs at large x is essential for determining 
high-energy cross sections at collider energies such as in 
search of new physics beyond the Standard Model, where 
structure information at large x feeds down through per- 
turbative Q 2 evolution to lower x and higher values of Q 2 
Q. Thus, it is of paramount importance to decrease the 
uncertainties in the extraction of PDFs by deriving the 
parameterizations directly from data without resorting 
to theoretical assumptions alone for the extrapolation to 
x ~ 1. 

Extending to larger x at a finite Q 2 means encounter- 
ing the resonance region. An important consequence of 
duality is that the resonance and deep inelastic regions 
are deeply connected and properly averaged resonance re- 
gion data could facilitate our understanding of the deep 
inelastic region. In some of the QCD analysis performed, 
the higher-twist terms have been extracted from data se- 
lected with a kinematic cut of W 2 > 10 GeV 2 jJOM HI • 
However, it was shown in several analyses 24, 25} that 
only a relatively small higher-twist contribution consis- 
tent with the one obtained in Ref. [3, [H, H3 would 
be necessary to describe the entire F 2 structure function 
spectrum. Indeed, S. Liuti et al. analyzed resonance 
region data within a fixed-VF 2 framework. This study 
found that the higher-twist contributions in the reso- 
nance region are similar to those from W 2 > 10 GeV 2 , 
with the exception of A region where the effects seem 
to be larger. This is in no way surprising if one thinks 
of it as a consequence of quark-hadron duality which en- 
sures that, on average, higher- twists are small or cancel. 
As a consequence of duality, the wealth of resonance re- 
gion data could be utilized to access the large x region 
and constrain the PDFs in this regime. This approach, 
however, requires a very good understanding of the Q 2 
dependence of the data in these kinematic regions of x 
and Q 2 where the perturbative evolution is no more the 
only mechanism responsible for the Q 2 behavior. 

Figures 1 and 2 depict three pQCD parameterizations 
of the F2 structure function at two Q 2 values, 4 and 
8 GeV 2 , with a zoom- in of the large-x region. The 
CTEQ6M parameterization shown is a QCD fit to hard 
scattering and DIS data (BCDMS, NMC, CCFR, E605, 
CDF, HI, ZEUS, DO) with Q 2 > 4 GeV 2 and W 2 > 12.25 
GeV 2 . The x dependence of the PDFs is parameterized 
at a Q 2 of 1.3 GeV 2 and then the QCD evolution equa- 
tions are utilized to evolve the distributions at higher Q 2 
in the NLO (and LO). The authors employed the twist- 
2 pQCD formalism so the kinematic cuts used for data 
selection were tested to ensure that the introduction of 
simple phenomenological higher twist terms would not 
improve significantly the quality of the fit. The CTEQ6M 
fit shown in Figs. 1 and 2 is obtained in the MS (mod- 
ified minimal subtraction) factorization scheme. One of 
the main improvements over earlier CTEQ fits is the ad- 
dition to the global set of data of new measurements 
(HI, ZEUS, DO) which provide better constraints on the 



PDFs, in general, and on the gluon distribution at large 
x, the result being a harder gluon distribution in this 
region. The other noteworthy improvement is the full 
treatment of uncertainties of the PDFs and their physi- 
cal predictions, using an eigenvector-basis approach. 

The MRST2004 parameterization is a QCD fit to a 
wide set of deep inelastic and related hard scattering 
data (BCDMS, SLAC, NMC, CCFR, CDF, HI, ZEUS, 
HERA, DO) with Q 2 > 2 GeV 2 and W 2 > 12.5 GeV 2 . 
The x dependence of the PDFs is parameterized at Q 2 
of 1 GeV 2 and a fixed order, LO or NLO or NNLO, 
QCD evolution is performed to specify the distributions 
at larger Q 2 where data exist. A global fit to the data 
then determines the parameters of the input distribu- 
tions. Though the fits are performed in the standard 
MS scheme, the gluon distribution is parameterized in 
the DIS (deep-inelastic scattering) factorization scheme 
and then transformed to the MS scheme. Together with 
more precise calculations of the splitting functions up to 
NNLO, this is actually the main improvement over earlier 
MRST fits (MRST2001 [H). Indeed, the NLO global 
analysis with this new gluon parameterization appears 
to work extremely well when compared to Tevatron jet 
data and is even better for the NNLO fit. This objective 
couldn't be accomplished by previous MRST parameter- 
izations. 

Both CTEQ6M and MRST2004 are shown here with 
target mass corrections (TMC) included according to 
Ref. [13]. CTEQ6M has more strength at large x than 
MRST2004. For most part, this discrepancy originates 
from the fact that the two groups use different functional 
forms for the parameterization of the non-perturbative 
input parton distributions and neither parameterization 
is constrained by measurements in the large x regime. 

The ALEKHIN parameterization shown in Figs. 1 and 
2 is an update of an earlier parameterization [14j |. the 
significant improvement being the use of recent calcu- 
lations of the exact NNLO evolution kernel. The data 
used were from SLAC, BCDMS, NMC, HERA, HI, ZEUS 
with kinematic cuts of Q 2 > 2.5 GeV 2 , W 2 > 3.24 GeV 2 
and x < 0.75. The model for the data description was 
based on pQCD with phenomenological parameterization 
of the twist-2 and higher-twist contributions to the struc- 
ture functions. The analysis was performed in the MS 
scheme with the number of flavors fixed at 3. The twist-2 
PDFs were parameterized at Q 2 = 9 GeV 2 . The pQCD 
analysis was done up to NNLO. Given the rather low W 2 
cut used to select the data set, ALEKHIN parameteriza- 
tion includes, besides the typical parameters of pQCD, 
parameters to account for the target mass and dynam- 
ical higher-twist effects. This is a novelty, considering 
the standard procedure of performing QCD fits. The 
higher-twist contributions to the structure function were 
parameterized in additive form: 



F 2 = F 2 



LT.TMC 



H 2 (x) 



(4) 



LT TA4C 

where F 2 ' has contributions from the twist-2 terms 



6 



with target mass corrections included according to Ref. 
(27| and the dynamical twist-4 term Hi (x) is parameter- 
ized in a model-independent way as a piece-linear func- 
tion of x. The use of the exact NNLO corrections made 
possible an improvement in the positivity of the gluon 
distributions extrapolated to small x and Q: in this pa- 
rameterization the gluon distributions are positive up to 
Q 2 = 1 GeV 2 , i.e. throughout the kinematic region where 
the parton model proved to be applicable. Since the 
ALEKHIN parameterization is based on fits to data with 
lower W 2 than CTEQ6M and MRST2004, its PDFs are 
expected to be better constrained at large x. 

Parameterizations of Ff. The parameterizations 
discussed above provide parton distribution functions 
from which the nucleon structure function can be con- 
structed in QCD frame-work. ALLM97 is a fit to just 
the nucleon (proton) structure function. In order to 
construct the structure function for a nucleon inside 
a nucleus substantial additional challenges need to be 
overcome. There are a host of well-documented is- 
sues in extracting nucleon structure functions from nu- 
clear data, even from deuterium data (see, for instance 
[H [H, Ha HH H2). At large x in particular, the effects 
of Fermi motion, nuclear binding, the EMC effect, off- 
shell corrections, and the like are quite large, and must 
be taken into account. Since there is no consensus on how 
best to accomplish this, we have here chosen to compare 
the measured deuterium resonance region data directly 
to deep inelastic deuterium structure functions. Specifi- 
cally, we have chosen to multiply the array of structure 
functions previously discussed by the following parame- 
terization of d/p (deuteron over proton) 33]: 

— = 0.9851 - 0.5648a; - 0.0904a; 2 + 0.7183a; 3 - 0.3428a; 4 . 
2p 

(5) 

This equation is the result of a data fit up to x — 0.8, and 
may not be constrained correctly at the highest x. More- 
over, it assumes no Q 2 -dependence, a need for which has 
been indicated in other works (3~ij . A more thorough ap- 
proach might in the future consider specifically structure 
functions formed from nuclear PDFs, such as those found 

in m,[3i[3l[3|. 

Target Mass Corrections. At large enough values 

of Q 2 and W 2 , QCD provides a rather clear and rigorous 
perturbative description of the physics that generates the 
Q 2 behavior of the structure function. When W — > M, 
where M is the proton mass, both the nonperturbative 
kinematical power corrections (target mass corrections) 
and the dynamical higher twist have to be taken into 
account. Since these characterize the long-range non- 
perturbative interactions between quarks and gluons, the 
dynamical higher-twist terms contain information about 
the dynamics of confinement. The target mass correc- 
tions, on the other hand, arise from purely kinematic ef- 
fects associated with finite values of Q 2 jv = AM 2 x 2 /Q 2 . 
The target mass terms are related to the twist-2 opera- 
tors and contain no additional information on the non- 
perturbative multi-parton correlations. In consequence, 



TABLE I: The kinematic regime covered by E00-116 at a 
beam energy of 5.5 GeV. 



Angle(deg 


) Momentum(lieV/c) 


x • 


Q ((^ eV ) ) 




2.26 






37.93 


1.94 
1.0/ 


0.48-0.92 


3.58-5.48 




2.17 






41 


1.86 
1.60 


0.53-0.94 


3.99-5.86 




1.94 






45 


1.67 
1.44 


0.55-0.95 


4.28-6.29 




1.34 






55 


1.16 
1.47 


0.60-0.94 


5.01-7.07 




1.31 






60 


1.19 
1.04 
0.89 


0.52-0.95 


4.52-7.38 


70 


0.91 
0.80 


0.60-0.83 


5.38-7.11 



target mass effects should either be corrected for in the 
data or the effect should be included in the QCD fits if 
one aims for a consistent comparison of data to QCD 
fits. The target mass effects were taken into account in 
the CTEQ6M, MRST2004 and ALEKHIN parameteriza- 
tions of the structure function according to the prescrip- 
tion of Georgi and Politzer [27|. It is non-trivial to note 
that there is not an univers ally agreed-upon prescription 
to account for target mass [39|, HO, HU , and so the choice 
of approach inherently introduces some uncertainty to 
this analysis. 



III. EXPERIMENT AND DATA ANALYSIS 

The experiment E00-116 was carried out in Summer 
2003 in Hall C, at JLab. A fixed electron beam of energy 
5.5 GeV came incident on cryogenic targets. The target 
system consisted of 4 cm long liquid hydrogen and deu- 
terium, contained in circular aluminum cans. Scattered 
electrons were detected in the High Momentum Spec- 
trometer (HMS). The Short Orbit Spectrometer (SOS) 
was used for detection of positrons, which was used to 
estimate possible electron background originating from 
charge-symmetric processes such as ir° production and 
subsequent decay in the target. The data were taken at 
various scattering angles and momenta as follows: for 
each fixed spectrometer angle, the central momentum 
was varied in order to cover a region in W 2 from about 1.2 
to 4.5 GeV 2 . The kinematics covered by this experiment 
are shown in Table 1. These data extend the existing 
Hall C resonance region measurements at larger x and 
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A. Experimental Setup 

1. Beam Line 

During E00-116, the Continuous Electron Beam Accel- 
erator Facility (CEBAF) at JLab provided an unpolar- 
ized, Continuous Wave (CW) electron beam of 5.5 GeV, 
with currents up to 100 /J.A. The beam was steered from 
the Beam Switch Yard to the experimental hall through 
the beam line. Hall C beam line is equipped with mag- 
nets used to focus and steer the beam, as well as several 
monitors needed to measure the energy, current, position 
and profile of the beam. The profile and the absolute po- 
sition of the beam is monitored utilizing superharps. A 
superharp consists of a frame and three tungsten wires 
(two horizontal and one vertical) which are moved back 
and forth through the beam to determine the centroid 
position to about 10 ^m. However, the superharp can- 
not be used during the data taking because it has a de- 
structive interaction with the beam. Therefore, Beam 
Position Monitors (BPMs) [42| are used to continuously 
monitor the relative beam position during data taking. 
The BPMs are nondestructive to the beam and are cal- 
ibrated with superharp scans. During this experiment, 
the typical relative variation of the beam position at the 
target was found to be less than 0.2 mm. 

The beam energy is measured using the superharps 
and the dipole magnets in the beam line. Due to the fact 
that the dipole fields are accurately mapped and that 
the beam path is determined with high precision by the 
superharps, the accuracy of the absolute beam energy 
measurement is at the level of 5 x 10~ 4 GeV. 

The beam current and charge in Hall C is measured 
by a system of beam current monitors (BCMs) together 
with a parametric current transformer (Unser) [431 ] . All 
these monitors are placed in the beam line before the 
target in the following order: BCM1, Unser, BCM2 and 
BCM3. Although the BCMs have a very stable offset, 
the gain drifts with time and the Unser is used for BCMs 
gain calibrations. Dedicated calibration runs are typi- 
cally performed to minimize the effects of drifts in the 
BCMs gains. For this experiment, BCM2 was used for 
monitoring but due to time constraints, no BCM cali- 
bration runs were taken. However, the experiment that 
ran just before E00-116 had the same set-point for the 
BCM2 gain, such that it was possible to use their cal- 
ibration runs taken 5 weeks and 1 week before this ex- 
periment |44|]. The difference in BCM2 gain and offset 
when using each run individually or the combined runs 
was below 0.1%, indicating that there were no significant 
drifts in the BCM2 gain over a month period. Figure 3 
shows that the difference between the current as given 
by the BCM2 after calibration and the current as given 
by the Unser (the residuals) are within 0.15% above 50 
/j,A. The calibration result from the combined runs was 
used to calculate the current and the charge for this ex- 
periment. For this experiment current regime, the nor- 
malization uncertainty in the current measurement was 
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FIG. 3: (Color online) The difference between the current as 
given by the BCM2 after calibration and the current as given 
by the Unser. The two calibration runs, 46203 and 46667 
were taken one month apart and yielded very similar results 
for the BCM2 gain and offset (see text). 

estimated to be ~0.3% at 100 /iA and originated from 
possible small drifts of the BCM gain, from the preci- 
sion of the BCM calibration and from the accuracy of 
the Unser in measuring the current, the latter bringing 
the largest contribution. The point-to-point uncertainty 
was estimated to be 0.05% by taking the difference in 
the normalization uncertainties propagated at 80 fj,A and 
100 /iA, this being roughly the range in the beam current 
used in this experiment. 

The electron beam generated by CEBAF is a high cur- 
rent beam with a very small transverse size (100-200 
FWHM). To prevent damage to the targets and to min- 
imize the changes in the cryotarget densities due to lo- 
calized boiling, a rastering system is used to distribute 
the deposited energy of the beam in an uniform man- 
ner over the target volume. The raster consists of two 
sets of steering magnets: the first set of magnets rasters 
the beam vertically and the second horizontally. For this 
experiment the raster consisted of a 2x2 mm uniform 
structure. A detailed description of the Hall C raster 
system is given in Ref . 45] . 

2. Target 

This experiment used the standard liquid hydrogen 
and liquid deuterium cryogenic target system in Hall 
C. The liquid targets were contained in aluminum cans. 
Data were taken on aluminum foils (dummy target) for 
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background measurement and subtraction. Also several 
runs were taken using a carbon foil target in order to 
determine the beam offsets relative to the pivot of the 
target. The cryogenic target cells were mounted on a 
cryostack together with the combination of carbon and 
aluminum target sled. The cryogenic system ensured 
that the temperature and density of the liquid targets 
were maintained during data taking at optimum values 
of 19 K and 0.0723 g/cm 3 for hydrogen and 22 K and 
0.1674 g/cm 3 for deuterium. 

To accurately determine the experimental luminosity 
it is necessary to have a precise knowledge, among oth- 
ers, of the targets' density and thickness. Though the 
cryogenic system is designed to ensure that the liquid hy- 
drogen and deuterium targets are maintained at a fixed 
nominal temperature in all conditions, in reality, when 
the beam passes through the targets and deposits heat 
there are local changes in the temperature and density of 
the cryogen (boiling effect). Dedicated data (luminosity 
scans) are taken to study and correct for this effect. Dur- 
ing E00-116, luminosity scans were performed on both 
hydrogen and deuterium targets. It was found that the 
boiling effect gives a small correction of (0.35 ± 0.32)% 
/ 100 fiA to the luminosity for both cryogenic targets 
as seen in Fig. 4. This parameterization was utilized 
on a run-by-run basis to correct for the boiling effect for 
both liquid targets. The majority of the data were col- 
lected at ~100 fiA. The density correction at this current 
is of the size of the uncertainty of the fit, therefore the 
normalization uncertainty was taken to be 0.35%. The 
point-to-point systematic uncertainty on the density cor- 
rection, originating from the uncertainty in the current, 
is negligible. 

Because of the circular geometry of the cryogenic tar- 
get cell, a careful analysis is required to determine the 
effective target length that will enter in the calculation 
of the luminosity. If the beam was exactly aligned along 
the diameter of the target, the effective target length will 
simply be the outer diameter of the target cell minus the 
cell walls. If there is a displacement between the beam 
and the center of the target, then the effective target 
length will be — dx 2 , where r is the inner radius 

of the cell and dx is the beam offset from the center of 
the target. For E00-116, several sources of information 
were used to determine the effective target length: the 
target survey that provides measurements, at room tem- 
perature, of the outer diameter of the cell together with 
the thickness of the cell walls, a survey of the target posi- 
tion relative to the pivot, and dedicated data taken on a 
central Carbon foil that provide the beam offset relative 
to the pivot The effective target length used in the 
cross section extractions is listed in Table 2. 



3. Spectrometers 

In what follows, a summary of the main characteris- 
tics of Hall C spectrometers will be given with emphasis 
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FIG. 4: (Color online) Relative hydrogen and deuterium tar- 
get yield versus beam current. The correction for the boiling 
effect is obtained from a linear fit to the data shown as the 
dotted line. 



TABLE II: E00-116 effective lengths of the cryogenic targets. 



Target Effective Target Length(cm) 



Hydrogen 
Deuterium 



3.946 ± 0.029 
3.927 ± 0.029 



on the aspects relevant to E00-116. Detailed information 
about the Hall C HMS and SOS can be found in Ref. 
[47l [48[ . The HMS is a magnetic spectrometer consisting 
of a 25 deg vertical bend dipole magnet for momentum 
dispersion and three quadrupole magnets for focusing. 
All magnets are superconducting. For this experiment, 
the HMS was operated in the point-to-point optical tunc. 
The range used in the momentum (E 1 ) acceptance, S = 
-2-, was of ±8% while the range in the angular (8) ac- 
ceptance, A(8) was ±35 mrad. The SOS consists of a 
quadrupole magnet and two dipole magnets. For E00- 
116, the point-to-point optical tune was used. The range 
used in S was of (-15, +20)% while the range in the an- 
gular acceptance was ±60 mrad. 

The detector packages for the two spectrometers are 
very similar and consist of two drift chambers for track re- 
construction, scintillators arrays for triggering, a thresh- 
old gas Cerenkov and an electromagnetic calorimeter, 
which were both used in this experiment for particle iden- 
tification (PID) and pion rejection. The HMS Cerenkov 
counter was used to distinguish between e~ and ir~ with 
momenta between 0.8 and 2.3 GeV. For this purpose, the 
Cerenkov tank was filled with Perfluorobutane (C4F10, n 
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Number of Photoelectrons 

FIG. 5: An example of the distribution of number of photo- 
electrons collected in the HMS Cerenkov detector. The pion 
peak appears at zero while electrons produce on average about 
13 photoelectrons. As it will be discussed later in the text, 
a cut of number of photoelectrons = 2 was used to separate 
electrons from pions. 



= 1.00143 at 1 atm and 300K) at about 0.9 atm making 
the detector fully sensitive to e~ but insensitive to tt~ in 
the momentum range specified above. The SOS Cerenkov 
counter was used to detect e + and reject ir + with mo- 
menta ranging from 0.8 to 1.7 GeV. The Cerenkov tank 
was filled with Freon-12 (n = 1.00108 at 1 atm) at about 

I atm giving a Cerenkov threshold of 3 GeV for n + and 

II MeV for e+. For E00-116, a typical spectrum of the 
HMS Cerenkov number of photoelectrons is shown in Fig. 
5. In this distribution, the ir~ events peak at zero while 
the e~ events give an average of about 13 photoelectrons. 

The HMS and SOS calorimeters are identical except 
for their size. Each calorimeter consists of 10 cm x 10 
cm x 70 cm blocks of TF-1-1000 type lead glass (p = 3.86 
g/cm 3 , n = 1.67 and radiation length r.l. = 2.5 cm). The 
HMS calorimeter is 13 blocks high while the one in the 
SOS just 11. The calorimeters are rotated by 5 deg from 
the spectrometers optical axis in order to reduce eventual 
losses through the cracks between the blocks. For this ex- 
periment, the calorimeters were used to detect e~ (HMS) 
and e + (SOS). The hadrons that could have reached the 
calorimeters were mostly 7r~ or 7r+. The e~ (e + ) were 
distinguished from 7r~ (n + ) according to their fractional 
energy, total energy deposited in the calorimeter, normal- 
ized by the momentum. The e~ (e + ) deposit their entire 
energy in the detector peaking in the fractional energy 
spectrum at 1; the 7r~ (tt + ) deposit around 0.3 GeV and 
they will peak in the fractional energy distribution at 0.3 
GeV/E . A typical distribution of the fractional energy 
deposited in the HMS calorimeter is shown in Fig. 6. 
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B. Data Analysis 

The inclusive electroproduction cross section can be 
expressed as: 



d 2 a 
dQ.dE' 



(N, 



measurei 



-BG) 



1 1 1 



N e N t dttdE' A e ' 



(6) 



Here, N e is the number of incident electrons and N t is 
the number of target particles per unit area, which can 
be calculated in terms of the mass density p, the atomic 
number A and the thickness x from N t — pN £ x (Na is 
Avogadro's number). N measurec i is the number of scat- 
tered electrons observed in the solid angle dQ, and in the 
energy range dE . BG is the background, A is the detec- 
tor acceptance and e is the detector efficiency. The most 
significant corrections that were applied to N measurec i will 
be discussed below. 



FIG. 6: An example of the distribution of fractional energy 
deposited in the HMS calorimeter. Electrons deposit their en- 
tire energy in the calorimeter peaking at 1 in the distribution 
while pions loose just a fraction of their energy. In the data 
analysis, a cut of fractional energy deposited = 0.7 was used 
to separate electrons from pions. 



1. PID Cut Efficiency 

The rejection of negatively charged pions was accom- 
plished by placing requirements on both the number of 
photoelectrons collected by the Cerenkov detector, num- 
ber of photoelectrons larger than 2, and the fractional 
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energy deposited by the particle in the calorimeter, frac- 
tional energy larger than 0.7. In what follows, the effi- 
ciency of these cuts in not rejecting valid electrons will 
be discussed. 

Cerenkov Cut Efficiency. In order to determine 
how many electrons are lost when applying the Cerenkov 
cut number of photoelectrons larger than 2, it is impor- 
tant to work with a clean sample of electrons (no pion 
contamination). Once a clean sample of electrons is se- 
lected, than the Cerenkov cut efficiency is determined 
from the ratio of number of events that pass the cut to the 
total number of events in the clean sample. If the sample 
is pion contaminated then the Cerenkov cut efficiency will 
be artificially lower. Unfortunately, it was impossible for 
this experiment to select a clean sample of electrons just 
with a calorimeter cut. For the particular kinematics of 
E00-116, the pion to electron ratio was rather large (up to 
150:1). The pions can undergo charge- exchange reactions 
and deposit up to their entire energy in the calorimeter. 
This could result in a high-energy tail for pions which 
could extend beyond 1 in the fractional energy spectrum 
making the selection of a clean sample of electrons prac- 
tically impossible even with high calorimeter cuts. The 
unbiased electron cut efficiency for the Cerenkov was de- 
termined by extrapolating to zero pion to electron ratio. 
It was found to be (99.60 ± 0.24)%, which in very good 
agreement with the findings of other experiments that 
ran at similar conditions [4{|. This value was used as 
a correction for the data. The normalization systematic 
uncertainty was taken to be the uncertainty of the fit 
extrapolation at zero pion to electron ratio, 0.24%. 

Calorimeter Cut Efficiency. Just as for the 
Cerenkov, the estimation of how many valid electrons are 
lost when using a cut on the fractional energy deposited 
in the calorimeter was complicated by the fact that, for 
this experiment, the large pion to electron ratio made im- 
possible the selection of a clean sample of electrons using 
just a cut on the number of photoelectrons acquired in 
the Cerenkov detector. The same approach was taken in 
this case as for the estimation of the Cerenkov cut ef- 
ficiency: the calorimeter cut efficiency was extrapolated 
to zero pion to electron ratio in order to find the true 
electron efficiency. The extrapolation was done for each 
momentum setting separately in order to deconvolute the 
efficiency dependence on the pion to electron ratio from 
the dependence on the resolution of the calorimeter. An 
example of the efficiency extrapolation at zero pion to 
electron ratio for one momentum setting is shown in Fig. 

7. The efficiency, obtained in this manner, and pa- 
rameterized as a function of momentum is shown in Fig. 

8. The parameterization was used as a correction in the 
data analysis. Our parameterization was compared to 
the result obtained from an experiment that ran in sim- 
ilar experimental conditions but at different kinematics 
where the pion to electron ratio is small [49']. The two 
results were found to agree within 0.3%. The normaliza- 
tion and point-to-point uncertainties on this correction 
were estimated to be 0.3% and 0.25%, respectively. 
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FIG. 7: (Color online) The electromagnetic calorimeter cut 
efficiency versus n/e ratio. The fit represented by the solid 
line extrapolates the efficiency at zero 7r/e ratio in order to 
obtain the true electron efficiency. 
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FIG. 8: (Color online) The electromagnetic calorimeter cut 
efficiency versus momentum. The solid line is a fit to the 
data and parameterizes the cut efficiency dependence on the 
momentum of the particle. This parameterization was used 
to correct for the loss of valid electrons due to the calorimeter 
cut inefficiency. 
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2. Backgrounds 



There are three physical processes that are possible 
sources of background for this experiment: electrons scat- 
tered from the target aluminum walls, negatively charged 
pions that are not rejected by the PID cuts, and electrons 
originating from other processes like charge symmetric 
processes which produce equal number of positrons. Each 
of these possible sources of background will be discussed 
in what follows. 

Target Cell Background. During data taking on the 
cryogenic targets, some of the incoming electrons scatter 
on the aluminum walls of the target cell and end up being 
detected at the same kinematics as the electrons that 
scatter from the cryogen. This background has to be 
determined and subtracted from the measured yields in 
order to obtain the yields for scattering from the cryogen 
only. To determine this background dedicated data were 
taken on a dummy target at exactly the same kinematics 
as on hydrogen and deuterium. In order to minimize the 
data acquisition time, the total thickness of the dummy 
target was about eight times the total cell wall thickness 
seen by the beam. The background coming from the 
scattering from the target cell walls BTW(E ,9), was 
determined as: 



BTW(E 



(7) 



where 9 is the spectrometer angle, Q w (d) is the to- 
tal charge incident on the cell walls (dummy), T w r^\ is 
the total thickness of the cell walls (dummy), Nd [E , 
0) is the number of events collected from the dummy 
run after applying efficiencies and dead time corrections 
and R^u) 1S t ne external radiative correction (external 
bremsstrahlung emission) for the cell walls (dummy). 
The target cell background subtraction was performed 
for each hydrogen and deuterium run on a (E ,9) bin- 
by-bin basis. The size of this background was at most 
18% and its uncertainty was dominated by the statistical 
uncertainty on Nd(E , 9) and by the uncertainty in mea- 
suring the thickness of the cell walls and dummy. The 
thickness of the cell walls was known up to 1% while, 
by comparison, the uncertainty in the dummy thickness 
measurement was negligible [501 . This led to a systematic 
uncertainty in the cross section of at most 0.2%. The sta- 
tistical uncertainty on Nd(E , 9) was propagated to the 
uncertainty of the cross section. 

Pion Background. Even after applying the PID 
cuts, some pion background may still be present. Al- 
though pions do not produce Cerenkov light directly, they 
can generate, through ionization, (5-rays in the materi- 
als preceding the Cerenkov detector (electron knockout). 
These knock-on electrons could have high enough energy 
to emit Cerenkov light and pass the PID Cerenkov cut. 
In the electromagnetic calorimeter, the pions give sig- 
nal according to their energy loss but through a charge- 
exchange reaction they can produce neutral pions which 
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FIG. 9: (Color online) The pion contamination (see text) as 
a function of momentum. As expected, the pion contamina- 
tion is larger for the deuterium target (bottom panel) than 
for the hydrogen target (top panel). The solid lines represent 
parameterizations of the pion contamination as a function of 
momentum. These parameterizations were used as correc- 
tions in the cross section extraction. 



decay into 77 or e~e + 7. In this way, the entire energy 
of the pion can be deposited in the calorimeter. This 
process typically gives the high-energy "tail" for pions 
that extends to deposited fractional energy of 1. For this 
experiment, the pion background was estimated using a 
method developed for the Hall C E99-118 analysis [34| 
in which the pion rejection factor is used to normalize 
the pion fractional energy distribution in the calorime- 
ter. The number of events in this normalized distribution 
that pass the PID calorimeter cut of 0.7 represents the 
pion contamination. The result of the pion contamina- 
tion estimation for this experiment is shown in Fig. 9 for 
both hydrogen and deuterium targets. As expected for 
a heavier target, the pion contamination for deuterium 
(maximum of about 1.7%) is about 3 times larger than 
for hydrogen (maximum of 0.5%) at the lowest momen- 
tum. For this background subtraction, parameterizations 
as a function of momentum were used as corrections (no 
angle dependence was observed) . The point-to-point un- 
certainty on this correction was determined to be 0.2% 
for both targets. 

Charge-Symmetric Background. For electron- 
proton scattering there is a significant probability to pro- 
duce a neutral pion in the target which then decays into 
77 or e~e + 7. The photons can further convert into 
electron-positron pairs in the target material or in the 
materials preceding the detectors. Photons can also be 
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produced through the Bethe-Heitler process. However, 
the leptons resulting from Bcthc-Heitler processes are 
very forward peaked so their contribution is significant 
only at forward angles. The outcome is that the sec- 
ondary electrons will end up being detected together with 
the scattered electrons. For the kinematics of this ex- 
periment (backward angles) the dominant source of sec- 
ondary electrons is the neutral pion production in the 
target and the subsequent decays. The electron pro- 
duction through Bethe-Heitler process is negligible. The 
fact that the background electrons are produced in pairs 
with positrons (charge-symmetric background) can be 
exploited experimentally and the background electrons 
can be disentangled from the scattered electrons by de- 
tecting positrons. 

For E00-116, due to the limited running period (less 
than a week), it was decided to take advantage of the 
availability of the SOS. The SOS has a larger momentum 
acceptance than the HMS and two SOS momentum set- 
tings could easily cover an HMS scan with three momen- 
tum settings. However, by using a different spectrom- 
eter for positron measurements the photon to electron- 
positron pairs conversion factor is different as the pho- 
tons encounter different radiation lengths of material. In 
addition, the SOS acceptance function and the detector 
inefficiencies are different than the HMS ones. Taking 
into account these considerations, we decided that an ac- 
curate estimation of the charge-symmetric background 
would require the extraction of the positron cross sections 
rather than the yields as it was done for previous Hall 
C experiments. This way, the charge-symmetric back- 
ground would be corrected by subtracting the measured 
positron cross section from the measured electron cross 
section bin-by-bin on a (E ,8) grid. 

The first step in the positron cross section analysis was 
to perform the detectors calibrations. Once the calibra- 
tions were performed, the positron yield selected with 
PID cuts was binned in the (6,6) acceptance around the 
central values. The yield was corrected for the electronic 
and computer dead times and for the tracking ineffi- 
ciency. In order to obtain the positron yields from the 
cryotargets alone, the endcap contributions had to be 
subtracted. Also the pion contamination was determined 
and parameterized as a function of momentum for each 
cryotarget and applied as correction to the yield. Next, 
the spectrometer acceptance corrections were calculated 
and applied to the yield. Thus, the positron cross section 
was obtained on a (6,8) grid. 

Our goal was to determine the cross section at the cen- 
tral angle as a function of momentum but still to keep 
the statistics accumulated. This could be done, in prin- 
ciple, by statistically averaging the measured cross sec- 
tion over the angular acceptance. However, the varia- 
tion of the positron cross section across 8 acceptance was 
non-negligible. Therefore, before averaging, a model was 
needed to remove the cross section 8 dependence (the so 
called bin-centering correction). The positron cross sec- 
tion model used for this purpose was developed by P. 
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FIG. 10: (Color online) An example of the measured positron 
cross section across the SOS angular acceptance. The cross 
section is shown at various stages in the analysis. In black 
circles the cross section is depicted before acceptance and 
bin-centering corrections. The blue squares show the cross 
section after acceptance corrections were applied. The angu- 
lar dependence of the positron cross section is obvious. The 
red triangles represent the cross section after both the accep- 
tance and bin-centering corrections were applied. It can be 
seen that the bin-centering corrections removed the angular 
dependence of the cross section. 



Bosted |5l|. The model uses a fit to the charged pion 
production data accumulated at SLAC [52| . The neutral 
pion production is estimated as the average of the pos- 
itive and negative pion production. The positron cross 
section is calculated using the decay branching ratios for 
a neutral pion and the radiation length of the material 
where a photon that results from the decay can produce 
electron-positron pairs. Taking into account that the 
positron cross section model described above was used 
just for bin-centering, the main requirement was that the 
model should describe the shape of the angular depen- 
dence of the positron cross section. To make sure that 
this requirement was met, first it was checked if, after ap- 
plying the bin-centering correction, there is any angular 
dependence left across the acceptance. A typical exam- 
ple is shown in Fig. 10 where it can be seen that, within 
the statistical uncertainty, the bin-centering correction 
removes the angular dependence of the cross section. 

It was also checked that the data overlap in the angular 
acceptance from one central angle setting to the next if 
the model would be used to center the data at certain an- 
gle values in the acceptance. Good agreement was found 
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FIG. 11: (Color online) The measured positron cross section 
on hydrogen as a function of momentum compared to the 
model of P. Bosted [H]]. 



FIG. 12: (Color online) The measured positron cross section 
on deuterium as a function of momentum compared to the 
model of P. Bosted [H. 



for neighboring scans in the overlapping region. 

Finally, the positron cross section was extracted at 
fixed central angles as a function of momentum. Fig- 
ures 11 and 12 show the positron cross section for both 
hydrogen and deuterium targets compared to the model 
of P. Bosted. It can be seen that the model describes 
qualitatively well the momentum dependence of the cross 
section. 

As stated previously, the electron data, both the scat- 
tered and the background electrons, were taken using 
HMS while the charge-symmetric background was mea- 
sured using SOS. Therefore, at the end of the experiment, 
a setting was taken at the same kinematics in HMS and 
SOS (both spectrometers were set on negative polarity). 
The result of the analysis of this scan in the two spec- 
trometers is shown in Fig. 13. It was found that the anal- 
ysis in the two spectrometers agreed within 1.3%. This 
translated in a normalization uncertainty in the scattered 
electron cross section below 0.2%, considering that the 
relative contribution of the charge symmetric background 
to the measured cross section was at most 15%. 

In the end, the charge symmetric background was sub- 
tracted bin- by-bin in (E ,8). For the subtraction, the 
positron cross section had to be centered at the scat- 
tered electron data kinematics. The model of P. Bosted 
was used for this purpose. Quantitatively, it was found 
that the use of the model for bin-centering corrections 
in momentum introduced an uncertainty of 6% in the 
positron cross section at 38 deg, 41 deg, 45 deg, 55 deg 
and of 20% at 60 deg and 70 deg. This translated in an 



uncertainty on the electron cross section up to 2% at the 
lowest momentum at 60 deg and 70 deg, but below 0.2% 
for the rest of the data. 



3. Acceptance Corrections 

For a fixed angle and momentum setting, the spec- 
trometers have a finite angle and momentum acceptance. 
This experiment used the same procedure of extracting 
the spectrometer's acceptance functions as previous Hall 
C experiments. This procedure is described in detail else- 
where [53j . 

For this analysis, the acceptance correction was ap- 
plied on a bin-by-bin basis in (6,8). The point-to-point 
uncertainty on the acceptance correction in HMS was es- 
timated to be 0.8%. This is dominated by the position 
uncertainties on the target, collimator, magnets, and de- 
tector package. The normalization uncertainty on the 
acceptance correction was determined by combining in 
quadrature an uncertainty of 0.7% coming from the re- 
duction in the solid angle and an uncertainty of 0.4% due 
to modeling of the HMS optics [53| . 



C. Cross Sections Extraction 

For this experiment, the electroproduction differential 
cross section H(e, e') and D(e, e') was extracted accord- 
ing to Eq. 6, binned in 16 and 20 bins in momentum 
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FIG. 13: (Color online) The comparison of SOS and HMS 
analyses for H(e,e') (top panel) and D(e,e') (bottom panel). 



and angle, respectively. Our goal was to obtain the one 
photon-exchange (Born) cross section at a fixed angle as 
a function of momentum. For this, two additional correc- 
tions were necessary: the bin-centering correction, which 
makes possible the extraction of the differential cross sec- 
tion at a fixed angle without sacrificing statistics and the 
radiative corrections which are necessary to obtain from 
the measured cross section the one photon-exchange con- 
tribution. Considering that both the bin-centering cor- 
rections and the radiative corrections are calculated using 
a model for the cross section, the sensitivity of our results 
to the model input was studied in detail. All this will be 
discussed in what follows. 

Bin-Centering Corrections. As previously men- 
tioned, the measured cross section was initially extracted 
binned in small momentum and angle bins corresponding 
to the acceptance intervals in S and 0, respectively. The 
goal, however, was to extract the cross section at a fixed 
angle, the central angle 9 C , keeping all the statistics ac- 
cumulated. If the cross section would not vary across the 
9 acceptance, then the cross section at 9 C could be sim- 
ply obtained by statistically averaging the cross section 
over the angular acceptance. However, for the kinemat- 
ics of this experiment, the variation of the cross section 
over the angular acceptance was not negligible. Thus, 
before statistically averaging, the so called bin-centering 
correction had to be applied in order to center the cross 
section measured in the 9 acceptance interval at 9 C . This 
correction was applied as: 



a data (E, E , 9 C ) = a data (E, E , 9,) 



data 



a model (E,E ,9 C ) 
a model {E,E' ,0iY 



(8) 

where <j model (E, E' , 9. t ) and a model (E, E' ,9 C ) are the 
model cross sections calculated at 9 L and C , while 
a data (E,E' ,9 t ) and a data (E, E' ,9 C ) are the cross sec- 
tions extracted from the data at 9i and 9 C . The bin- 
centering correction was applied to the measured radi- 
ated hydrogen and deuterium cross sections using the ra- 
diated model cross section. The models used to calculate 
the correction will be discussed next. 

Radiative Corrections. In the perturbative pic- 
ture, the lowest order process in a (the electromagnetic 
running coupling constant) that contributes to the cross 
section for inclusive electron-nucleon scattering is repre- 
sented schematically in Fig. 14 (a). Besides this leading 
one photon exchange diagram (Born), there are higher 
order processes in a that contribute to the scattering. 
These diagrams are shown schematically in Fig. 14 (b, c, 
d, e) and include vacuum polarization (the exchanged 
photon creates particle-antiparticle pairs), vertex pro- 
cesses (emission and reabsorption of virtual photons), 
and Bremsstrahlung (emission of real photons in the field 
of the nucleon during interaction) . In order to determine 
the differential cross section that accounts just for the 
one photon exchange process, all the other contributions 
from higher order processes in a have to be calculated 
and corrected for in the measured cross section. The ra- 
diative processes can be divided into two main categories: 
internal and external. The internal effects take place at 
the scattering vertex and include Bremsstrahlung, vac- 
uum polarization, vertex processes and multiple photon 
exchange. External Bremsstrahlung occurs within the 
target material before or after the primary scattering 
takes place and is dependent on the target thickness. As 
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FIG. 15: Schematic representation of the processes that can 
contribute to the E00-116 measurements. 

a consequence, the energy of the incoming and/or the 
scattered electron will change. 

The measured cross section can be written as the sum 
of various processes (Fig. 15): 

hydrogen radiated ■ radiated /n\ 

measured inelastic ' elastic ' \ J 

deuterium radiated _r_ radiated , radiated (~\(\\ 

measured inelastic ' elastic ' quasielastic \ J 

To obtain the Born inelastic cross section, the radiative 
tails from elastic/nuclear elastic and quasielastic cross 
sections were subtracted while the inelastic radiative ef- 
fects were corrected multiplicatively. In practice, the 
usual calculation of radiative corrections includes only 
the emission of one hard photon. However, there is a 
probability for the electron to emit two hard photons. 
Therefore additional corrections (the a 2 term) should be 
taken into account when estimating the radiative correc- 
tion. 

There are two programs utilized for the radiative cor- 
rections calculations: one based on the Mo and Tsai for- 
malism I54| and the other one based on the Bardin for- 
malism [55j . which includes in the calculations the two 
hard-photon radiation and has a different treatment of 
the soft photon contribution. The program based on Mo 
and Tsai formalism calculates both internal and exter- 
nal radiative corrections, unlike the Bardin one which 
calculates just internal radiative corrections. For this ex- 
periment, the radiative corrections (both internal and ex- 
ternal) were estimated using the Mo and Tsai formalism. 
The a 2 term correction was estimated using the Bardin 
formalism and at this experiment kinematics it proved to 
be below 0.5% (1%) for hydrogen (deuterium). The size 
of the correction is within the theoretical uncertainty of 
the calculation and it was assigned as a point-to-point 
uncertainty. 

Additional uncertainties in the experimental cross sec- 
tion originate from the elastic (elastic and quasielastic) 
contribution subtractions. At the kinematics of this ex- 
periment, the clastic contribution to the total experimen- 
tal cross section is negligible. The uncertainty coming 
from the quasielastic contribution subtraction was esti- 
mated by propagating the point-to-point model uncer- 
tainty into the experimental cross section. This kine- 
matic dependent uncertainty was parametrized for each 
angle setting separately. Finally, the normalization un- 
certainty in the cross section coming from the theoretical 



uncertainty in the radiative corrections calculation was 
estimated to be 1% (56j . 

Iteration Procedure As stated previously, the bin- 
centering corrections were calculated using a model for 
the cross section. The same model was typically used to 
calculate the radiative corrections. In order to minimize 
the model dependence of the extracted cross section, an 
iterative procedure was followed. First, a starting model 
was used to calculate both the bin-centering and radia- 
tive corrections. Then, the extracted cross section was fit 
and the new fit was used to calculate the corrections and 
reextract the cross section. This process continued until 
the extracted cross section did not vary significantly (not 
more than 0.3%) from one iteration to the next. Addi- 
tionally, the iteration procedure was followed using two 
different starting models. After the last iteration, the 
two sets of cross sections were expected to be consistent. 

For H (e, e') the two starting models (fits to previous 
data) used were the model of M.E. Christy and P. Bosted 
[5^ and H 2 model [HH . The fitting procedure used is de- 
scribed extensively in Ref. (5?| • Only two iterations were 
necessary and it was found that the difference in the cross 
section between the last iteration and the one before last 
was about 0.3%. Also after each iteration the difference 
in the cross section when starting with the two mod- 
els specified above was calculated. This difference after 
the last iteration was assigned as a kinematic dependent 
uncertainty accounting for the model dependence of the 
final result. 

For D(e,e') measurements, the iteration was per- 
formed using the same fitting procedure as for H(e, e') 
but with one modification: the fit form for the resonances 
used non-relativistic Breit-Wigners. The data seemed to 
be described better around pion threshold by such a fit. 
Three iterations were performed and the difference in the 
cross sections between the last iteration and the second 
last was around 0.3%. Just like for the H (e, e') data set, 
two different models were used in the iteration proce- 
dure: the Bodek model [59[ and ALLM97 [To| multiplied 
by the parameterization of the ratio of the deuterium and 
the proton electroproduction cross sections [H[ . The dif- 
ference between the two sets of cross sections after the 
last iteration was parameterized to give the uncertainty 
originating in the possible model dependence of the final 
result. 

Systematic Uncertainties The total point-to-point 
systematic uncertainty in the cross section extraction was 
taken as the sum in quadrature of the individual un- 
certainties. An overview of these uncertainties is given 
in Table 3. The total normalization systematic uncer- 
tainty amounted to about 1.75% while the statistical un- 
certainty is below 3%. The total and Born deradiated 
differential cross sections extracted from this experiment 
for both H(e,e') and D(e, e'), together with the point- 
to-point associated uncertainties are given in Tables 5-9. 
The extracted Born differential cross sections are also 
shown as a function of x and W 2 in Figs. 16-19. These 
data provide large x and intermediate Q 2 high-precision 
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FIG. 16: (Color online) The H(e,e') Born differential cross 
section extracted from E00-116 data at a beam energy of 5.5 
GeV and spectrometer central angle of 38 deg (empty cir- 
cles), 41 deg (empty squares) and 45 deg (empty triangles) 
as a function of W 2 (top panel) and x (bottom panel). Both 
the statistical and point-to-point systematic uncertainties are 
plotted. The curves shown represent the fit after the last 
iteration 46, H3. 



FIG. 17: (Color online) The H(e,e') Born differential cross 
section extracted from E00-116 data at a beam energy of 5.5 
GeV and spectrometer central angle of 55 deg (empty cir- 
cles), 60 deg (empty squares) and 70 deg (empty triangles) 
as a function of W 2 (top panel) and x (bottom panel). Both 
the statistical and point-to-point systematic uncertainties are 
plotted. The curves shown represent the fit after the last 
iteration [46, WA- 



TABLE III: Point-to-point systematic uncertainties in the ex- 
perimental parameters and the corresponding systematic un- 
certainties in the differential cross section. 



Quantity 


Uncertainty 


S a {%) 


Beam Energy 


5xl0 -4 


0.30% 


Scattered e Energy 


5xl0 -4 


0.25% 


Scattered e Angle 


0.2 mrad 


0.26% 


Beam Charge 


0.05% 


0.05% 


Dead time 


0.25% 


0.25% 


Trigger Efficiency 


0.2% 


0.2% 


Tracking Efficiency 


0.2% 


0.2% 


PID cut efficiency 


0.25% 


0.25% 


Pion Cont. Subtraction 


0.2% 


0.2% 


Charge- Symmetric Background 


6% - 20% 


<2% 


Acceptance Correction 


0.8% 


0.8% 


Radiative Corrections 


0.5%-3.6% 


0.5%-3.6% 


Model dependence 


0.2%-5% 


0.2%-5% 



measurements in the resonance region where the preci- 
sion of existing data from SLAC is typically 5 to 30% for 
the statistical uncertainty alone. This precision is not 
enough to distinguish between theoretical parameteriza- 
tions of the structure function which at a Q 2 of 8 GeV 2 
and x = 0.8, for example, differ by at most 30% as shown 
in Fig. 2. 



D. F 2 Extraction 

The structure function F 2 was calculated utilizing the 
formula: 



F 2 = 



d 2 a 1 + R Kv 1 



1 



dVLdE' 1 + eR Ana T l 



(11) 



where K = (W 2 - M 2 )/(2M), v = E - E and a is the 
electromagnetic coupling constant. The quantity e is the 
degree of polarization of the virtual photon: 



-tan 2 -)- 1 , 



(1 



and r is the flux of the virtual photon: 



aK E 1 



2ir 2 Q 2 El-E 



(12) 



(13) 



The extraction of F 2 requires the knowledge of both the 
differential cross section and the quantity R which is the 
ratio of the longitudinal to the transverse component of 
the cross section. For this experiment, it was not planned 
to measure R but just the differential cross section since 
R is not expected to be large in the Q 2 range of this data 
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FIG. 18: (Color online) The D(e, e') Born differential cross 
section extracted from E00-116 data at a beam energy of 5.5 
GeV and spectrometer central angle of 38 deg (empty cir- 
cles), 41 deg (empty squares) and 45 deg (empty triangles) 
as a function of W 2 (top panel) and x (bottom panel). Both 
the statistical and point-to-point systematic uncertainties are 
plotted. The curves shown represent the fit after the last 
iteration [4rj |. 



set. Thus, for the structure function extraction, it was 
needed to resort to an existing R parameterization. In or- 
der to estimate the sensitivity of the Fi structure function 
to different R parameterization, the following approach 
was taken. F2 was extracted at the lowest and highest 
Q 2 for this experiment using three different R parame- 
terizations: R from R1990 [gfl], R1998 [H[ and E94-110 
(56| . respectively. It was found that F2 varies on aver- 
age by 2% when different R parameterizations are used. 
R1990 and R1998 are parameterizations of R extracted 
mainly from DIS measurements, while R from E94-110 
is extracted from resonance region measurements, only. 
The R parameterization from E94-110 is kinematically 
limited to W 2 < 3.85 GeV 2 and typically lower Q 2 than 
this experiment. However, it was shown to agree, where 
applicable kinematically, to R1990 and R1998. The R 
parameterization R1998 was obtained using a larger data 
set than R1990 (R1998 used, in addition to the data set 
of R1990, other measurements extending the parameter- 
ization to lower and higher x) , and also had a better con- 
fidence level of the fit (73%) than R1990 (61%). For all 
these reasons, R from R1998 was used for the structure 
function calculation and an additional uncertainty of 2% 
was assigned to F2 in order to account for the sensitivity 
of the extraction to R. 




FIG. 19: (Color online) The D(e, e') Born differential cross 
section extracted from E00-116 data at a beam energy of 5.5 
GeV and spectrometer central angle of 55 deg (empty cir- 
cles), 60 deg (empty squares) and 70 deg (empty triangles) 
as a function of W 2 (top panel) and x (bottom panel). Both 
the statistical and point-to-point systematic uncertainties are 
plotted. The curves shown represent the fit after the last 
iteration [4^]. 

IV. RESULTS 

In this section we present our results with a preamble 
of previous quark-hadron duality studies. With this very 
brief summary of previous studies we intend to create 
the apropriate context for the detailed discusssion of our 
results and the conclusions that will be drawn. 



A. Previous Quark-Hadron Duality Studies with 
Electron Scattering 

Over three decades ago, Bloom and Gilman ac- 
knowledged the resonance-scaling connection in inclusive 
electron-nucleon scattering. After more than 20 years, it 
was an early JLab experiment that revived the interest 
in the phenomenon of quark-hadron duality @, [62|, [63| • 
This experiment confirmed the observations of Bloom 
and Gilman and, in addition, acknowledged the onset 
of duality also locally. The findings of these studies 
prompted interest in a more detailed analysis of duality, 
one within the QCD formalism. 

Comparisons of the resonance region data to some 
QCD predictions were performed using measurements 
from JLab experiment E94-110 @, [Hj]. The F i and F l 
extracted from E94-110 were compared to QCD fits from 
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the MRST [H[ and CTEQ6 [H| collaborations evaluated 
at the same Q 2 as the data and with the inclusion of tar- 
get mass corrections. It was observed that the QCD fits 
seem to describe on average the resonance strength at 
each Q 2 value investigated - Q 2 of 1.5, 2.5 and 3.5 GeV 2 - 
as if the resonances would follow, on average, the same 
perturbative Q 2 evolution as the QCD fits. On a more 
quantitative level, quark-hadron duality was investigated 
by computing ratios of the integrals of the structure func- 
tion over x in the resonance region at fixed Q 2 values from 
both the data and the scaling curves. It was found that 
duality seems to hold better than 5% above Q 2 — 1 GeV 2 
when compared to MRST with target mass corrections. 
However, at the highest Q 2 of 3.5 GeV 2 and the highest 
x the ratio to MRST was noticed to rise above unity up 
to 18% 0. 

This finding came in strong contradiction with the ex- 
pectation that duality should work even better at higher 
Q 2 . If the higher-twist contributions seem to be small or 
cancel to some degree at low Q 2 then, considering that 
these terms are weighted by powers of 1/Q 2 in the op- 
erator product expansion, it is expected that this must 
be even more the case at a higher Q 2 . Moreover, the 
observation of increasing discrepancy between data and 
some QCD fits with increasing Q 2 (and increasing x) is 
not unique to the resonance region: DIS data from SLAC 
exhibit the same behavior [67| . In consequence, this rise 
has been ascribed not to a violation of duality but rather 
to an underestimation of the large x strength in some 
QCD fits. 

These studies made obvious the utility of high preci- 
sion resonance region data at an even higher Q 2 (and 
thus larger x). This extension of resonance region mea- 
surements was crucial for the verification of QCD fit be- 
havior in this kinematic regime. Considering that most of 
the currently available large x data lie in the resonance 
region, the confirmation of quark-hadron duality as an 
effective tool would offer much needed experimental con- 
straints for theoretical predictions in the region of x — > 
1. 

In what follows, therefore, we present quark-hadron 
duality studies performed using the F 2 structure function 
extracted from this experiment, as well as from earlier 
Hall C [5l,[62] and SLAC [H M, M, [El III M, III 
measurements in the resonance region. 

B. Quark-Hadron Duality: the Q 2 Dependence 

An exhaustive description of nucleon structure in terms 
of parton distribution functions requires knowledge of the 
strength of the PDFs for the entire x regime. Most global 
QCD fits are essentially unconstrained at large x [3^ |. 
The quark-hadron duality phenomenon could be the key 
for providing experimental constraints in the large x re- 
gion by the use of properly averaged resonance data. This 
avenue relies, however, on our ability to unravel the Q 2 
dependence of the data in a region where the perturba- 



tive mechanisms are not the only ones to be taken into 
account. 

In this context, a comparison of the Q 2 dependence 
of various theoretical predictions to the one exhibited by 
averaged resonance re gion data was studied in a similar 
fashion to that of [56|, l62| . Ratios of the integrals of the 
F 2 structure function were considered: 

F% aram (x,Q 2 )dx' [ ' 

The integrand in the numerator -P& ata - is the F 2 struc- 
ture function extracted from the experimental cross sec- 
tions. The integrand in the denominator _pP aram -_ [ s the 
F 2 structure function as given by the parameterizations 
introduced in Sect. 2: CTEQ6M+TM, MRST2004+TM, 
ALEKHIN, ALLM97. It is important to note that, for 
this analysis, F. 'P aram - was generated at the same values 
of x and Q 2 as the data and was integrated over the 
same range in x as the data using the same integration 
procedure. This, by dint of the W 2 cuts used to obtain 
these global fits, by definition extends them into regions 
where they are not constrained in x, and only their Q 2 
dependence is determined. 

For global duality studies, the limits of the integrals, 
x m i n and x max , were the experimental x values corre- 
sponding to W mm = 1.3 GeV 2 and W^ ax = 4.5 GeV 2 , 
respectively. To compare the Q 2 dependence of theo- 
retical predictions to individual resonance structures, for 
local duality studies, the resonance regions were delim- 
ited using the same W 2 cuts as in a previous analysis 
@: 

1. first region (I s *) — > W 2 G [1.3 , 1.9] GeV 2 

2. second region (2 nd ) — > IV 2 € [1.9 , 2.5] GeV 2 

3. third region (3 rd ) —> W 2 £ [2.5 , 3.1] GeV 2 

4. fourth region (4 th ) -> IV 2 e [3.1 , 3.9] GeV 2 

5. DIS region (DIS) -> IV 2 e [3.9 , 4.5] GeV 2 

These W 2 limits translate in the integrals of Eq. 14, to 
x m i n and Xmax values according to: 

X= W 2 + Q 2 - M ' (15) 

where M is the proton mass. As an example, the x range 
covered by different resonance regions for two Q 2 values 
are given in Table 4. At a given Q 2 , the lowest W 2 re- 
gion (the first region) corresponds to the highest x regime 
while for a fixed W 2 region, the larger the Q 2 , the larger 
the x regime. 

Figures 20-23 depict the results of the global and local 
duality studies performed for H(e,e'). The quantity / 
is shown for each resonance region individually as well 
as integrated over the full region specified above. The 
uncertainties shown are obtained by adding in quadra- 
ture the statistical and systematic uncertainties on the 
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TABLE IV: An example of the x ranges covered by different 
resonance regions for two Q 2 values. 





x range 


W z region 


Q A = 2 (GeV 2 ) 


Q' J = 6 (GeV") 




0.66 - 0.83 


0.85 - 0.93 


cyTbCt 


0.55 - 0.66 


0.79 - 0.85 




0.47 - 0.55 


0.73 - 0.79 




0.40 - 0.47 


0.67 - 0.73 


DIS 


0.35 - 0.40 


0.62 - 0.67 



numerator alone. No parameterization uncertainties are 
plotted. The latter are typically substantial at large x, 
on the order of 100% [ll[. Overall, the data of this exper- 
iment (blue circles), previous JLab data (the data rep- 
resented by black stars and red triangles are from [62| 
and 15611 . respectively), and SLAC data (green squares 
[H, M, H Hi, M, III III III), are found to be in 
good agreement. A slight disagreement could be ob- 
served between this experiment and the SLAC experi- 
ment E-8920 which is singled out from the other SLAC 
data sets (empty green square). The disagreement be- 
comes smaller as we approach the DIS region. This is 
possibly related to the fact, that for E-8920, the radia- 
tive corrections were rescaled to bring the DIS data in 
agreement with the other SLAC experiments which could 
however have resulted in a possible incorrect estimation 
of the radiative corrections in the resonance region. 

Figure 20 presents our results when we compare the 
integral of the F^ extracted from the data to the integral 
of F 2 P obtained from the CTEQ6M PDFs with the inclu- 
sion of target mass effects as explained in Sect. 2. The 
quantity / is close to unity at a Q 2 of about 1.5 GeV 2 , 
and then rises above unity with increasing Q 2 . However, 
/ reaches a plateau at a Q 2 of about 4 GeV 2 and, above 
this value, the Q 2 dependence saturates. This behavior 
is displayed when the integration is done globally as seen 
in the bottom right panel, but also for "all" the individ- 
ual resonance regions except for the first resonance region 
(upper left panel) . The saturation of the Q 2 dependence 
indicates that the discrepancy between data and param- 
eterization is not a Q 2 dependent effect. It is most likely, 
therefore, due to the fact that CTEQ6M+TM does not 
model accurately the strength of the PDFs at large x. 
Put differently, / being greater than 1 and, reaching a 
constant value above Q 2 ~ 4 GeV 2 most likely does not 
represent the failure of QCD in describing the Q 2 evo- 
lution of the averaged resonance region data but rather 
a paucity in the strength of the PDFs at large x. The 
resonance region data do display on average a QCD type 
Q 2 -dependence. 

The ratio / seems to become constant at a slightly dif- 
ferent value for each resonance region. In fact, as we 
move from the fourth resonance region (/ ~ 1.1) to the 
third (/ ~ 1.28) and then to the second (/ ~ 1.35) the dis- 
crepancy increases. This is possibly related to the grow- 
ing uncertainty associated with PDF strengths at large 
x. For a fixed Q 2 , the second resonance region probes a 



larger x regime than the third and, the third larger than 
the fourth, etc. Since the x dependence of the PDFs is 
less and less constrained at larger and larger x, we expect 
this to be reflected in a more obvious way when we study 
the second resonance region than the third, for example. 

It was reported before that the AT- A transition region 
provides a different behavior when compared to the rest 
of the resonance region [24| ■ This could be related to the 
fact that this region is the only one with a single resonant 
state and there are arguments that more than one state 
is necessary to approximate closure and duality [76j]. It 
should also be pointed out that the first resonance region 
probes the highest x regime where the PDFs are expected 
to be least constrained. 

Figure 21 shows the ratio of the integrals of F% from 
the data and MRST2004 with target mass corrections. 
The observed Q 2 dependence of I yields similar conclu- 
sions to those drawn from the comparison to CTEQ6M: 
we encounter the same rise of I with Q 2 which eventually 
saturates for all resonance regions except for the first one 
and also globally. And just as for CTEQ6M, I saturates 
at a different values for each resonance region. This is 
not surprising considering that the extraction procedure 
of PDFs for MRST2004 is rather similar to the one em- 
ployed by CTEQ6M. There are, however, few features 
that set apart the comparison with MRST2004. This pa- 
rameterization undershoots the data by an even larger 
amount and I saturates at a larger value of Q 2 than for 
CTEQ6M+TM. This most likely results from the differ- 
ences in the x dependence modelling of the PDFs between 
the two parameterizations. 

Figure 22 shows the comparison of our averaged F% 
resonance data to the parameterization of ALEKHIN. 
It should be reminded that it is not just the leading 
twist that is considered in this parameterization, as is 
the case with CTEQ6M and MRST2004, but also the 
higher- twists. By explicitly accounting for higher- twist 
terms, Alekhin can extend the validity of his fit to an x as 
large as 0.75 and a W 2 as low as 3.24 GeV 2 . Though this 
W 2 cut practically excludes resonance region data, is still 
more permissive than the cuts employed by CTEQ6M or 
MRST2004, ensuring that the ALEKHIN fit is far better 
constrained at large x. Indeed, the agreement between 
the averaged resonance data and ALEKHIN is obvious 
in Fig. 22. For the fourth resonance and the DIS regions 
(upper and middle right panel, respectively) the quantity 
I is very close to unity across the entire Q 2 range inves- 
tigated. Good agreement is obtained when analyzing the 
second and third resonance regions: / deviates from unity 
by only 5% or less and, for most part, seems independent 
of Q 2 . This finding is quite remarkable: according to 
ALEKHIN higher-twist coefficients in the resonance re- 
gion, on average, differ from the ones extracted from the 
DIS region by at most 5%. It should be pointed out that 
not all of this already small discrepancy can be attributed 
to the contribution from higher- twist terms: ALEKHIN 
does not include resonance data in his fit therefore the 
x-dependence of the PDFs is unconstrained in this re- 
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gion, though to a far lesser extent than for CTEQ6M 
or MRST2004. This finding is quite different from what 
was observed when comparing to CTEQ6M, for exam- 
ple. There / deviates from unity by about 10% in the 
DIS region but by almost 35% in the second resonance 
region. When compared to ALEKHIN, the first reso- 
nance region (upper left panel) behaves differently but in 
this kinematic regime the fit validity is questionable. The 
data are well described on average also globally. Thus the 
higher-twist terms contributions in the resonance region 
is shown to be quantitatively comparable on average with 
the ones extracted from the deep inelastic scattering data 
pointing to the onset of quark-hadron duality. 

The comparison of our integrated F% resonance data to 
ALLM97 is presented in Fig. 23. The Q 2 dependence of I 
shows very good agreement between data and this param- 
eterization in the fourth resonance (upper right panel) 
and DIS (middle right panel) regions and also globally. 
If quark-hadron duality holds, this is to be expected con- 
sidering that ALLM97 successfully fits data down to a 
W 2 as low as 3 GeV 2 . The agreement slightly worsens 
as we move to the third and second resonance regions. 
Though / is about 7% above unity, it seems to be inde- 
pendent of Q 2 for the third resonance region. A familiar 
pattern emerges when analyzing the Q 2 dependence of / 
in the second resonance region: / rises with increasing Q 2 
but this rise eventually saturates around Q 2 of 4 GeV 2 . 
It is to be expected that as the larger x and lower W 2 
region is probed, the comparison between averaged reso- 
nance data and this parameterization to unravel some of 
its shortcomings like unconstrained x and Q 2 dependence 
or inability to fully account for target mass effects. The 
averaged resonance data in the first region (upper left 
panel) compare surprisingly well with ALLM97 but no 
definite conclusions could be drawn within QCD frame- 
work considering that ALLM97 accounts for the x and 
Q 2 dependence empirically and this is a region far from 
the domain of validity of this parameterization. Overall, 
the comparison of averaged resonance data to ALLM97 
confirms that, quantitatively, higher- twist terms contri- 
butions in the resonance region seem to be comparable, 
on average, with the ones in the deep inelastic scattering 
data. 

A similar pattern as for H (e, e') is observed when 
studying the Q 2 dependence of I for D(e,e') (Figs. 24- 
27). It should be reminded that, as discussed in Sect. 
2, there is an additional factor to consider when analyz- 
ing the results from our global and local quark-hadron 
duality studies for D(e,e'): all of the three QCD-based 
parameterizations utilized in our analysis provide PDFs 
from which the proton structure function is constructed. 
Also ALLM97 is a fit to only proton data. So, in order 
to obtain parameterizations for the deuteron structure 
function we used the F% parameterizations and the d/p 
parameterization from 33] introduced in Sect. 2. 

Figure 24 shows our results when we compare the in- 
tegral of F2 extracted from the data to the integral of 
F% obtained from CTEQ6M PDFs as explained above. 



The Q 2 dependence of I displays similar characteristics 
to the ones acknowledged in our study of the proton data 
presented in Fig. 20: / rises above unity with increasing 
Q 2 but a plateau is reached at a Q 2 of about 4 GeV 2 . 
Above this value / is practically independent of Q 2 . As 
observed before, the first resonance region (upper left 
panel) stands out and, in addition to the aspects dis- 
cussed for H(e,e'), the complication of having to resort 
to extrapolations of the d/p parameterization should be 
taken into account. 

A similar behavior is acknowledged when the data are 
compared with MRST2004 (Fig. 25). However, the Q 2 
dependence of I saturates at a larger Q 2 value than for 
CTEQ6M just as it happened for H(e,e'). This trend 
is even more accentuated as we probe larger x regimes 
(second resonance region in the middle left panel, for ex- 
ample) where the reliability of the d/p parameterization 
is questionable. 

Just as for 7J(e,e'), good agreement is observed when 
the D(e, e') data are compared to ALEKHIN, as seen in 
Fig. 26. In fact, except for the first resonance region 
(upper left panel) where both the PDFs and the d/p pa- 
rameterization are expected to be largely unconstrained, 
the -D(e, e') data are described by this parameterization 
down to the lowest Q 2 analyzed. Similar conclusions can 
be drawn from the comparison of the data to ALLM97 
which is presented in Fig. 27. 

To summarize, our studies showed that above a Q 2 
of about 4 GeV 2 for CTEQ6M and slightly higher for 
MRST2004 the ratio of the integrals of resonance data 
and parameterizations becomes independent of Q 2 . This 
is a very important finding which suggests that, above a 
surprisingly low Q 2 value, most of the disagreement be- 
tween the averaged resonance data and the above men- 
tioned parameterizations is unrelated to the violation of 
the Q 2 evolution by contributions from the higher-twist 
terms in the resonance region. In fact, the compari- 
son of our data to ALEKHIN and ALLM97 confirmed 
that higher-twist contributions to deep inelastic scatter- 
ing and averaged resonance region data are comparable. 
All these findings point to the unconstrained strength of 
the CTEQ6M and MRST2004 PDFs at large x as major 
source for the disagreement between data and the above 
mentioned parameterizations in this kinematic regime. 



C. Quark-Hadron Duality: the ^-Dependence 

Our quark-hadron duality studies discussed above in- 
dicate that there are small rather than large violations 
of the Q 2 evolution in the resonance region on average. 
Thus when referring to disagreements between data and 
theory, the ability of PDF-based calculations to describe 
the x-dependence of the data in particular at large x is 
brought into discussion. 

We used the averaged proton structure function data 
for the five W 2 regions to draw a comparison to the the- 
oretical calculations at fixed Q\ as a function of x. The 
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data averaging was done as follows: 

rx max F p dx 

-i-ip.ave JXmin 2, data tA a \ 

f p, = _ , (16) 

where x m i n and x rnax are the integration limits corre- 
sponding to the W 2 limits defined in the previous sec- 
tion. 

The averaged structure function data, F^ ave were then 
centered at a fixed Q\. Data within small Q 2 intervals 
were chosen for centering: for example, all F '^ ave data in 
the Q 2 interval of 2 GeV 2 to 4 GeV 2 were evaluated at 
Qq = 3 GeV 2 , etc. This was done as follows: 

Ql) = Fr ve (x, Q 2 ) F lvar am .^ Ql) 

The parameterization of M.E. Christy [l7| was used for 
bin centering because it describes the Q 2 dependence of 
the data to better than 3%. In addition, to study the 
sensitivity of the results to the choice of parameteriza- 
tion, the CTEQ6M fit was also utilized and the difference 
in the results when the two parameterizations are used 
was assigned as a systematic uncertainty. Just as for the 
quark-hadron duality studies discussed in the previous 
section our data were compared with all of the four pa- 
rameterizations introduced in Sect. 2 and the results are 
presented in Figs. 28-31. 

Given that the deviations between our locally av- 
eraged resonance region data and the expectations 
based on PDF parameterizations such as CTEQ6M and 
MRST2004 seem related to the uncertainty of these 
PDFs at large x, we start a comparison of the IE- 
dependence of our averaged resonance data with the phe- 
nomenological ALEKHIN and ALLM97 structure func- 
tion parameterizations. Figure 28 shows the ratio of the 
F% structure function extracted from the data as ex- 
plained above and the parameterization of ALEKHIN at 
four values of Q 2 as a function of x. At Q 2 of 3 and 5 
GeV 2 (upper right and lower left panel, respectively) the 
parameterization describes the x dependence of the data 
well, except for the largest x regime where measurements 
from the first resonance region are used. There is a small 
shift between the DIS+fourth resonance region data and 
the rest but no obvious disagreement which depends on 
x is observed. At Q 2 of 7 GeV 2 (lower right panel) our 
data probe the largest x regime where ALEKHIN is least 
constrained and we acknowledge a growing discrepancy 
between data and parameterization with increasing x. 
At Q 2 of 1 GeV 2 ALEKHIN fails to describe the IE- 
dependence of our data as x increases: the data probe 
here a regime where both the x and Q 2 limits of appli- 
cability are reached for this parameterization. Similar 
conclusions can be drawn from the comparison of our 
data with ALLM97 presented in Fig. 29. 

Figure 30 shows the comparison of the data to 
CTEQ6M+TM. The CTEQ6M uncertainties are plotted 
also as a band. The parameterization fails to describe 
the x-dependence of the data. The discrepancy is much 



larger than observed in our comparisons with ALEKHIN 
and ALLM97 and it also grows strongly with increasing 
x. The same conclusion can be drawn from the compar- 
ison of the data to MRST2004+TM presented in Fig. 
31. Again, this is not surprising considering that for 
both calculations the strength of the PDFs is largely 
unconstrained in the kinematic regime that we study. 
On the other hand, the ALEKHIN parameterization ac- 
counts explicitly for higher- twist, which allows it to in- 
clude data with lower W 2 (and larger x) than CTEQ6M 
and MRST2004. This offers better constraints to the x 
functional form which reflects in a more realistic descrip- 
tion of the data. 

Figure 32 shows the comparison of the F% structure 
function from Bourrely et al. [l5] to CTEQ6M+TM and 
ALLM97 at a fixed Q 2 value of 5 GeV 2 . The target 
mass effects where included in the parameterization of 
Bourrely et al. in an identical fashion as for CTEQ6M, 
according to Georgi and Politzer prescription (27J. At 
low x, x < 0.4, the three parameterizations agree reason- 
ably well. At large x however, significant discrepancies 
arise: up to 30% when compared to CTEQ6M+TM and 
larger for ALLM97, as shown in the insert of Fig. 32. 
Without the inclusion of target mass effects, the ratio 
of CTEQ6M+TM and ALLM97 to the parameterization 
from Bourrely et al. would be even larger at large x (x 
> 0.5). The major cause of this discrepancy is, most 
likely, the scarcity of high W 2 and high x data which 
could constrain the x dependence of the parton distribu- 
tion functions. The parameterization from Bourrely et 
al. undershoots the F^ structure function from ALLM97 
which was found to be in fairly good agreement with our 
averaged resonance region data. 

We conclude that what appeared to be a violation of 
quark-hadron duality when we compared our averaged 
resonance data with CTEQ6M and MRST2004 is actu- 
ally, for most part, a reflection of the inability of these pa- 
rameterizations to realistically model the large x strength 
of their PDFs. 



V. CONCLUSIONS 

We have performed high precision measurements of 
the H (e, e') and D(e, e') cross sections in the resonance 
region at large x and intermediate Q 2 . In this work, 
both global and local quark-hadron duality was quanti- 
fied for the proton and deuteron using our new large x 
data as well as previous resonance region measurements 
from JLab and SLAC. Previous studies @, HH indicated 
quark-hadron duality in the F^ structure function to hold 
better than 5% above Q 2 — 1 GeV 2 when compared to a 
typical QCD fit like MRST, with a growing discrepancy 
observed as regions of higher Q 2 , about 3.5 GeV 2 , and 
higher x were explored. This finding came in strong con- 
tradiction with the expectation that duality should work 
best with increasing Q 2 . The question arose whether this 
growing discrepancy was really a violation of duality by 
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contributions from higher-twist terms or mostly a conse- 
quence of the well-known issue of highly unconstrained 
PDFs at large x. 

We found that, when compared to CTEQ6M and 
MRST2004, the ratio of integrals of resonance data 
and these parameterizations becomes independent of Q 2 
starting with a value of about 4 GeV 2 for CTEQ6M, 
and slightly higher for MRST2004. This is an indication 
that, as expected from quark- hadron duality, there are 
only small violations of the Q 2 evolution by data in the 
resonance region, on average. This ratio saturates above 
unity at increasing values for regions with decreasing W 2 
(and increasing x), but remains constant in Q 2 , likely due 
to the uncertainty in the PDFs extraction at large x. 

The comparison to ALEKHIN revealed that the 
higher-twist contributions to the averaged resonance re- 
gion data are comparable to the ones in the low W 2 DIS 
region at the level of 5% or less. This points as well to the 
unconstrained PDFs at large i as a major source of the 
observed discrepancy between data and CTEQ6M and 
MRST2004. This argument is further supported by our 
studies of the x dependence of the data and parameteri- 
zations. 



This analysis concludes that, with a careful study of 
the Q 2 -dependent contributions, properly averaged reso- 
nance region data could be used to provide much needed 
constraints for PDFs at large x, shedding light on the 
parton dynamics in this regime. In view of quark-hadron 
duality, a CTEQ subgroup has begun to attempt the im- 
provement of PDFs at large x expanding the possible 
data sets by lowering the W 2 cut [73| . 
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FIG. 20: (Color online) Local and global duality studies for H(e,e') where CTEQ6M+TM [ill ] was used for comparison. 
Together with the data of this experiment (blue circles) are also plotted the results using measurements from two previous Hall 
C experiments, I. Niculescu (black stars) [§1 and E94-110 (red triangles) [56|, and from SLAC [H, [H 0, EL E2, El HI IH • 
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FIG. 21: (Color online) Local and global duality studies for H(e,e') where MRST2004+TM [ll| was used for comparison. 
Together with the data of this experiment (blue circles) are also plotted the results using measurements from two previous Hall 
C experiments, I. Niculescu (black stars) [g2 and E94-110 (red triangles) jH, and from SLAC [gl [H 0, El, El El 111, El • 
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FIG. 22: (Color online) Local and global duality studies for H(e,e') where ALEKHIN QJ, Q was used for comparison. 
Together with the data of this experiment (blue circles) are also plotted the results using measurements from two previous Hall 
C experiments, I. Niculescu (black stars) [§1 and E94-110 (red triangles) 0, and from SLAC [H [H El ED, El El El El • 
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FIG. 23: (Color online) Local and global duality studies for H (e, e') where ALLM97 [T(j was used for comparison. Together 
with the data of this experiment (blue circles) are also plotted the results using measurements from two previous Hall C 
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FIG. 24: (Color online) Local and global duality studies for D(e,e') where CTEQ6M+TM [l| multiplied by d/p ratio from 
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FIG. 25: (Color online) Local and global duality studies for D(e,e') where MRST2004+TM 12] multiplied by d/p ratio from 
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FIG. 27: (Color online) Local and global duality studies for D(e, e') where ALLM97 [T(| multiplied by d/p ratio from (3^] was 
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FIG. 29: (Color online) The ratio of F 2 structure function from data to F 2 from ALLM97 versus x at fixed Q 2 . 
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FIG. 30: (Color online) The ratio of F 2 structure function from data to F 2 from CTEQ6M+TM [Tl[ versus x at fixed Q 2 
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FIG. 31: (Color online) The ratio of F 2 structure function from data to F 2 from MRST2004+TM [T2| versus x at fixed Q 2 . 
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TABLE V: Differential cross sections extracted from the measurements of E00-116. The normalization uncertainty is 1.75%. 
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TABLE VI: Differential cross sections extracted from the measurements of E00-116. The normalization uncertainty is 1.75%. 
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TABLE VII: Differential cross sections extracted from the measurements of E00-116. The normalization uncertainty is 1.75%. 
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TABLE VIII: Differential cross sections extracted from the measurements of E00-II6. The normalization uncertainty is 1.75%. 
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TABLE IX: Differential cross sections extracted from the measurements of E00-116. The normalization uncertainty is 1.75%. 
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